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AN  ALGORITHMIC  SOLUTION  FOR  A  QUEUEING  MODEL  OF  A 


COMPUTER  SYSTEM  WITH  INTERACTIVE  AND  BATCH  JOBS 


ABSTRACT 

A  queueing  model  with  two  customer  classes,  one  with  infinite 
and  the  other  with  finite  source,  is  used  as  a  model  for  a  computer 
system  with  interactive  and  batch  jobs.  Using  an  imbedded  Markov 
Chain  representation  of  this  queueing  system,  and  an  algorithmic  ap¬ 
proach,  the  steady  state  joint  probability  distribution  of  the  number 
of  interactive  and  batch  jobs  at  a  job  service  completion  epoch  is 
derived.  Server  utilization,  mean  waiting  times,  joint  probability 
distribution,  and  mean  number  of  interactive  and  batch  jobs  at  an 
arbitrary  time  epoch  are  derived  using  these  probabilities,  discrete 
state  level  crossing  analysis  and  Little's  result. 


Key  Words:  Computer  Systems  Modelling,  Priority  Queues,  Imbedded 
Markov  Chain,  Algorithmic  Solution,  Level  Crossing 
Analysis. 
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INTRODUCTION 


The  intent  of  this  paper  is  to  present  an  algorithmic  solution 
for  a  queueing  model  of  a  computer  system  with  interactive  and  batch 
jobs.  Interactive  job3  arrive  from  a  finite  number  of  interactive 
terminals  and  batch  jobs  arrive  from  an  infinite  size  population. 

Once  an  interactive  job  arrives  to  the  computer  system,  the  corres¬ 
ponding  terminal  stays  passive  until  the  job  is  processed,  at  which 
time  the  terminal  becomes  active  and  begins  the  process  of  submitting 
a  new  job. 

The  queueing  models  associated  with  this  system  are  special  cases 
of  priority  queues  with  two  classes  of  customers  with  one  finite  and 
the  other  infinite  sources.  Models  in  which  the  finite  source  jobs 
having  preemptive  priority  over  infinite  source  jobs  have  been  exten¬ 
sively  analyzed  by  Avi-Itzhak  and  Naor  [1],  Colard  and  Latouche  [2], 
Jaiswal  [A],  Jaiswal  and  Thiruvengadam  [5],  and  Thiruver»gadam  [12]. 
While  the  technique  used  in  [1,  A,  5,  12]  is  complex  and  obtaining 
numerical  results  is  difficult,  the  algorithmic  approach  discussed 
by  Colard  and  Latouche  [2]  is  very  efficient.  This  efficiency  comes 
from  the  fact  that  the  above  model  can  be  treated  as  a  special  case 
of  an  M/PH/1  queue.  Simple  explicit  results  for  this  M/PH/l  queue 
can  be  found  in  Neuts  [10]. 

In  this  paper,  however,  we  shall  assume  that  the  priority  disci¬ 
pline  is  non-preemptive.  So  we  shall  model  the  computer  system  by  a 
queueing  system  with  two  classe  of  customers  and  non-preemptive  pri- 
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oricy  service  discipline.  Using  an  imbedded  Markov  chain  represen¬ 
tation  of  this  queueing  system  and  an  algorithmic  approach  discussed 
by  Lucantoni  and  Neuts  [8]  and  Neuts  [9],  the  steady  state  joint 
probability  distribution  of  the  number  of  interactive  and  batch  jobs 
at  a  job  service  completion  epoch  is  derived.  Server  utilization, 
mean  waiting  times,  joint  probability  distribution  and  mean  number 
of  interactive  and  batch  jobs  at  an  arbitrary  time  epoch  are  derived 
using  these  probabilities,  discrete  state  level  crossing  analysis  [11], 
and  Little's  result  [6]. 


2.  THE  MODEL 


B  Jfe; 


The  model  class  considered  here  is  of  a  single  server  queueing 
system  to  which  the  arrivals  form  two  independent  arrival  streams, 
one  from  a  finite  and  the  other  from  an  infinite  population  source. 
The  time  needed  by  each  of  the  M  members  in  the  finite  source,  to 
submit  a  job  (hereafter  called  interactive  job)  is  exponential  with 
mean  1/A^.  The  jobs  from  the  infinite  source  (hereafter  called  batch 
jobs)  arrive  to  the  server  according  to  a  Poisson  process  with  rate 
Aq.  The  service  requirements  of  these  jobs  at  the  server  are  general 
with  probability  distribution  function  Bq(*)  for  batch  and  B^(*)  for 
interactive  jobs  with  means  1/p^  and  1/p^,  respectively.  The  service 
discipline  is  non-preemptive  priority.  In  this  paper  we  consider 
two  cases,  one  in  which  the  interactive  jobs  have  higher  priority 
over  the  batch  jobs  and  the  other  in  which  these  priorities  are  re¬ 
versed. 


It  is  routinely  verified  that  the  number  of  interactive  and 
batch  jobs  juot  after  a  service  completion,  in  this  model,  forms  a 
Markov  chain.  Let  {(i,j),  i>0,  M>j^0>  represent  the  state  space 
with  i  representing  tb^  number  of  batch  and  j  representing  the  num¬ 
ber  of  interactive  jobs.  The  transition  probability  matrix  (TPM)  P, 
of  this  Markov  chain  is  of  the  form: 


where  the  matrices  A ,  r>0  and  r>0  are  all  square  substochasti< 


matrices  of  order  M+l.  The  matrices 


B  *  Z  £ 
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j . aartft . . . . . . 


are  both  stochastic. 


Let  Af(i,j)  and  B^i.j),  i  =  0,1,2,. ...M;  j  =  0,1,2, ... ,M,  be 
the  (i+l,j+l)th  entry  of  the  and  matrices,  respectively.  Then 
if  j^(n)  is  the  state  of  the  Markov  chain  at  the  n-th  transition, 

Ar(i,j)  -  P(£(n+1)  -  (k+r-l,j)|£(n)  =  (k,i)},  k  >  0 


and 


Br(i,j)  =  P{#(n+1)  =  (r,j)|$(n)  =  (0,i)} 


for  all  i,j  *=  0,1,2,. . .  ,M  and  n>0.  The  exact  values  of  these  entries 
will  depend  on  A^,  A^,  Bq(*),  B^(*)»  and  the  priority  assignment.  Be¬ 
fore  discussing  these  entries  any  further,  we  will  first  outline  the 
algorithmic  approach  used  to  analyze  the  Markov  chain  with  the  struc¬ 
ture  of  P  for  the  TPM.  This  approach,  as  a  general  methodology,  was 
introduced  by  Lucantoni  and  Neuts  [8]  and  Neuts  [9].  In  this  section, 
however,  we  will  present  the  necessary  results  only. 

Let  ^  be  the  steady  state  probability  vector  with  the  partion 
£  =»  (;£q,  ^,...)  where  i>0  are  all  (M+l) -vectors.  These  steady 
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'I 
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where  ^  »  (1,1,...)'  is  of  appropriate  dimension.  Next  we  present 
the  condition  for  the  existence  of  a  steady  state  solution.  Assum¬ 
ing  that  the  matrix  A  is  irreducible,  let  £  be  the  solution  to 


where 


and 


4<0)  •  rk 


4(n+1)  ‘  fek  +  E,4(n+1)AWr  +  Wn)AV 

r=2 


(6) 


fe’k  =  ($0Bk  +  «lAk)(I“Al) 


-1 


-1 


A'  =  A  (I-A.)  ,  r>0  . 

r  r  1  — 

These  iterations  of  x^Xn)  are  continued  until  a  sufficient  accuracy 
is  achieved  for  The  marginal  probabilities  are 


and 


(7) 


ci  =  ^ 


*  ’ 


for  batch  and  interactive  jobs,  respectively. 

It  can  be  shown  that  (see  equation  (40)  of  [8]) 


i =  s*1)  =  i*o  T-  Br">:1A'1(I"A+rt)"1  +  a~^h  (8) 

r=l 
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where  CZ)  =  £  jc  Zr  and  fl  is  an  (M+l)x(M+l)  matrix  with  identical 

r=l 

columns  equal  to  tt. 

The  first  moment  of  the  stationary  distribution  is 


&  « 


4.  ,■  T- 


w 


02 


JC(1)(1)  -  {~X(1)  [I-  £  rA]  +  I  Br 
r-1  r  ^  r-1  r 


+  ^0  S  rBr  •  >(I-A+n)_1  +  (£(1)  CD*)*. 

ral 


(9) 


where  ^C(l)  and  n  are  as  <  ifined  in  (8)  and  jg^(l)^  is  defined  in  the 
appendix. 

Next,  using  these  results  and  the  discrete  state  level  crossing 
analysis  [11],  we  will  derive  the  steady  state  joint  and  marginal 
probability  distributions  at  an  arbitrary  time  epoch. 

Let  g  be  the  steady  state  joint  probability  vector  at  an  arbitrary 
time  epoch.  Here  g  «  (g^g^...)  and  g±  “  (q.. j.qi;L» •  •  •  »qiM> •  Then, 
the  rate  of  upcrossings  from  the  compound  state  { (i, j ) ,  i<m,  j<n)  is 
(see  [11]), 


+  nA 


m 

1  l-o''1”' 


(10) 


since  the  arrivals  are  Markovian.  Now  let  x  ^  be  the  steady  state 
probability  that  a  type  r  (r*0  for  batch  and  r»l  for  interactive  job) 
departing  job  sees  i  batch  and  j  interactive  jobs  in  the  system. 

Then  the  rate  of  downcrossings  into  the  compound  state  ((i,j),  i<m, 
j<n)  is  (see  [11]), 


n  o  ®  1 

PnPn  r  +  piu-  1  x4~  » 
0  0  j-0  mj  1  A  i-0  in 


(11) 


where  p ^  is  the  fraction  of  time  the  server  is  busy  with  type  i  jobs. 
Now  equating  the  rates  of  up-  and  downcrossings,  from  (10)  and  (11) 
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we  get 


m 


m 


A E  a  .  +  nA,  I  q  =  p,yA  Ex  .  +  p,p,  E  x  . 

0  j-t)®3  1  i=0  in  0  0  j=0  mj  11  1=1  in 


Rearranging  terms,  we  get  a  recursive  formula  for  the  unknowns 

q_  .  It  is 
hnn 


%m  Ag+nA^P0p0x  mn  +  P1P1X  mn  +  ^g^p0p0x  mj 


m-1 


(12) 


+  ±*0  (plV  in  *  ">A,>  -  °in2>>  "E° 


Next  we  derive  the  marginal  pr./  -abilities  associated  at  an  arbitrary 
time  epoch. 

(1)  Interactive  Jobs 

Equating  the  rates  of  up-  and  downcrossings  (see  page  21  of  [11]), 
we  get 


nA-q1  =  p.p-x1  ,  0<n<M 

in  I  1  n  — 


(13) 


where  x  n  *  E  x  and  q  is  the  steady  state  probability  that  there 
i*0 

are  n  interactive  jobs  in  the  system  at  an  arbitrary  time  epoch.  Then 


1  plpl  1 
q  =  xx  ,  0<n<M 

n  n\,  n  ’  — 


3 

,5 


*  §?s 

# 

-5 
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(ii)  Batch  Jobs 

Using  a  similar  derivation  as  above  (see  [11]) ,  we  get 


i  0 

A_q 
0  m 


P0P0X  m  * 


(15) 


xm-  ^  x  n j  an<^  9  a  is  the  steady  state  probability  that  there 
are  m  batch  jobs  in  the  system  at  an  arbitrary  time  epoch.  Since  all 
batch  jobs  are  serviced,  PqUq  “  Aq.  Substituting  this  in  (15),  we 
have 


m 


x  ,  m>0 

m  — 


(16) 


The  following  section  provides  equations  to  evaluate  the  mean  number 
of  jobs  in  the  system  at  an  arbitrary  time  epoch,  and  the  mean  wait¬ 
ing  times. 


(i)  Interactive  Jobs 

Let  the  effective  interactive  job  arrival  rate  be  A1  .  Then 

e 

A*e  -  plwl*  **  is  the  “e®0  number  of  interactive  jobs  at  the 
server  at  an  arbitrary  time  epoch,  and  is  its  mean  waiting  time, 

\1e  -  (M-L1)A1 
or 

4  -  H  -  PjPj/Aj^  (17) 

Now  from  (17)  and  Little’s  result  [6], 
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Wx  =  (M/p1m1)  “  1/AX  . 


(18) 


(ii)  Batch  Jobs 

Let  LQ  and  be  the  mean  number  of  batch  jobs  in  the  system  at 
an  arbitrary  time  epoch  and  just  after  a  batch  jobs  service  comple¬ 
tion,  and  Wq  its  mean  waiting  time.  Then 


LQ  =  N°  ,  (19) 

since  =*  ,  and  from  Little's  result, 

m  m 

“o  ■  V\>  <20> 

To  use  (12),  (14),  (16)-(19),  and  (20),  we  need  xr^ ,  pr,  r  =  0,1, 
and  N^.  These  can  be  related  to  £  and  (1)  depending  on  the 
priority  rule  used.  We  shall  derive  these  relationships  later  for 
each  case,  separately.  The  following  definitions  will  be  used  in 
both  cases. 


|n,r:k)  **  P{n  batch  jobs  and  r  interactive  jobs  arrive  during  a 
time  interval  which  has  a  probability  distribution 
B^(*)»  given  there  are  k  interactive  jobs  left  in  the 

source  at  the  beginning  of  this  time  interval). 


10 


CASE  1:  Interactive  Jobs  Have  Higher  Priority  Over  Batch  Jobs 

Here  we  assume  that  the  interactive  jobs  have  higher  priority 
over  batch  jobs.  That  is,  if  there  is  at  least  one  interactive  job 
available  at  a  service  completion  epoch,  it  will  go  in  for  service. 
Then  the  entries  of  the  matrix  are  given  by, 

for  r>l 


f  P1(r-1, j+l-i:M-i)  M>i>0;  M>j>i-1 

Ar(i,j)  -  <  0  M>i>0;  i-l>j>0 

^  PQ(r,j:M)  i-0;  *£j>0 

and  for  r-0 

M>i>0;  *bj>0 

i-0;  M>j>0 


Similarly,  for  r>0 
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M>i>C;  M>j>0 


a0PQ(r, j:M)+a1P1(r,j :M-1)  ,  i«0:M-l>j>0 


^  aQP0(r,M:H) 


i®0;  j-M  , 


where  a, 


0  Xq+MX1 


xo  MXi 

and  a,  = 


1  xq+mx1 


(21) 


These  probabilities  can  be  evaluated  in  a  recursive  manner  and  the 
details  are  given  in  the  appendix.  Now  to  use  equations  (12) ,  (14) , 
(16)-(19),  and  (20),  we  need  jc1 ,  pQ,  p^,  and  N°.  We  will  derive 
these  next. 

Since  for  an  interactive  job  to  leave  behind  i  batch  and  n  inter¬ 
active  jobs,  an  interactive  job  must  go  in  for  service,  thus  the  use 
of  conditional  probability  and  convolution  arguments  lead  to 


i  n+1 
{  E  E 
k«l  j=l 


’V/i+l-k0  ,n) 


n+1 
+  E 

j-1 


x0jBi(.1,n)  +  a1xQOBi(0,n)}/p1 


and 


E  x 
i»0 


in 


(23) 


where  p^  is  the  probability  that  a  departing  job  is  interactive.  Now 
for  batch  jobs. 


(22) 
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!•  i 


0  ®+1 

x  mj  +  a0x00Bm^0,^^p0  * 


where  p  is  the  probability  that  a  departing  job  is  a  batch  job. 
Then 


0  “  0 
x  -Ex. 

“  j-0  mj 


(24) 


m+1  M 


M 


i  (.E.*10  .LVl-i<°>3)  +  a0x0Q  ,EA(0-J))  • 


■ o  i=*l  j-0 


j-0 


After  some  simplification,  and  using  the  normalizing  condition,  we  get 


nrt-1 

X"m  *  {  E  XiOP(nH'1~i)  +  a0X00P(°)}/p0  * 
i-1 


(25a) 


where  P(n)  is  the  probability  thf  n  batch  jobs  arrive  during  the 
service  of  a  batch  job,  and 


po  *  E  *10  +  Voo  • 

i-*l 


(25b) 


Using  the  theory  of  Markov  Renewal  Processes  (see  J.J.  Hunter  [3]), 
it  can  be  shown  that 


<i-lXi0  +  a°Xoo)  ,0/EC  * 


*  1  *  M  i  *oo 

EC  —  (  Ix..  +  anxnn)  7!  +  (  E  E  x. .  +  *iXnn)~—  + 


“OW  u0  \:0  J^ij  T  -roo^1  T  x0-mx1 
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and 


pl“  ( 


«  H 
Z  Z  x 
i-0  j-1 


+  Voo>^/EC- 


since  ECxEN  is  the  mean  busy  cycle  time  and  (  Z  x,n  +  a.xnr.)xEN/un  and 
•  H  i-1  0 

(  Z  Z  x. ,  +  a-xnn)xSN/w.  are  the  mean  time  in  a  busy  cycle  the  server 

i-0  j-1  J  1  UU  1 

is  busy  with  batch  and  interactive  jobs,  respectively,  where  EN  is  the 
mean  number  of  jobs  served  during  a  busy  cycle.  Note  that 


and 


.E  xi0  “  y0  “  X0C 
i-1 


«  M 

Z  Z  x  -  1  -  y 
i-0  j-1  13  U 


and  so  from  (25b) 


p0  *  y0  “  X00  +  a0x00 


(25c) 


Next  we  evaluate  N  .  Since 


N 


ZV 

m  x 

m-0 


m 


from  (25a)  we  get 
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00 


N 


0 


m+1 

E  m{  E  xi0P(nrfl-i)  +  x-^PCo)  }/pQ 
m»0  i*l 


00  00  CO 

-  {  E  Em  xi0P(m-(i-l))  +  xQ0a0  E  m  P(m)}/p0 
i*l  m*i-l  m*0 


*  {  E  ((i-1)  +  ,)  )x 
i-1 


iO 


+  xooVo}/p 


0 


CO  00 

{.5:ii  xiO  "  (1"pO^.E1xiO*  +  a0x00p0}</p0 
i=l  i**l 

(1)(0)  -  (1-pg)  ^o~X00^  +  a0x00p0^p0  * 


So 


L0  =  ^  “  (1“po^yo“xOO*  +  aoxoopo^po 


and  hence  *  L^/A^  can  be  obtained  once  we  evaluate  We  can  use 
(5)  and  (6)  to  obtain  this  provided  A  is  irreducible  and  inequality 
(A)  holds.  From  the  definition  of  ,  it  is  easily  verified  that  the 
(i+l,j+l)  the  entry  A(i,j)  of  A  is  given  by 


C  „  m_4  -A,t  ...  ,  -(M-j-l)A.t 

I  of  (j+l-i)(1"e  >j  e  dB1(t),  m>i>0,  m>j>i-l 


A(i,j)  -< 


oo  h  ~^1C  i  c 

0f  q>(l-e  1  ye 


1** 

A  dB0(t)  ,  i-0;  Mxj>0 


Otherwise 
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Similarly,  £  is  given  by 


(VvVvVv  ••••W 


Then  the  necessary  and  sufficient  condition  for  the  existance  of  ^  is 


p  =  ^  , 


where 


ttA 

•v 


If. 

'Ki 


CASE  2:  Batch  Jobs  Have  Higher  Priority  Over  Interactive  Jobs 

In  this  section  we  assume  that  the  batch  jobs  have  higher  pri¬ 
ority  over  interactive  jobs.  That  is,  if  there  is  at  least  one  batch 
job  available  at  a  service  completion  epoch,  it  will  go  in  for  ser¬ 
vice.  From  the  definition  of  P^(n,r:k)  it  is  easily  verified  that, 

for  r>0 


Ar(i,j) 


P0(r,j-i:M-i)  ,  M>i>0;  M>j>i 


M>i>0:  i>j>0 
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M>i>0 ; 


*0P0(r,j:M)  +  ^(r.jiM-l)  ,  i=0;M-l>J>0 


Br(i,j)  - 


a0P0(r,M:M) 


i»0;  j»M 


otherwise 


As  Indicated  in  the  previous  case,  these  probabilities  can  be  evalu- 

0  1 

ated  recursively.  Now  we  will  derive  expressions  for  Pq»P^»  ,  £  , 
and  Vp . 


Similar  to  the  previous  case,  using  Markov  Renewal  Theory,  we 


where 


°o  *  <  £  +  Voot*  > 

m=l  u 


^  1  00  i  *Q0 

EC  '  +  Vco^  +  +  VWuJJ +  x-ffijq 


pi  "  +  Voo)^/EC  • 


Note  that  r  (x  )  -  1  -  X0  and  E  x  -  x„  -  x  . 
m=l  j*l  J 


_ i 


The  relationships  for  %P  and  with  respect  to  ^  is  presented  next. 
These  results  follow  from  simple  convolution  arguments  and  condition¬ 


al  probabilities. 


m+1  j 


x  .  -  {  Z  l  x,.  Pn(mfl-i»  j-k:M-k) 
i«l  k=0  1K  U 


+  x00a0P0(m,j:M)}/p0  , 


where  p^  is  the  probability  that  a  departing  job  is  .  batch  job, 
and  so 


o  v  0 

m  j«0  “J 


i  2  x^  P(nrfl-i)  +  xooao  ?(“)  > /p0  » 
i*! 


where  p^  *  1  -  Xq  +  3qxqq.  Therefore 


»°  -  Z  m  x° 


{  Z  ra  Z  x^P(nri-l-i)  +  Z  a^^nPCm)  }/p. 
m«0  i*l  m*0  0 


{  Z  Z  mxiP(m-(i-l))  +  a0x00p0}/p0 
1=1  m=i-l 


{  z  (i-i-fpo)xi  +  Voopo}/po  • 
1*1 


(28) 


XD 


Lo  =  **  ^  po^1-xo^  +  aoxoopo^pc 


and  Wq  *  L^/Xq  can  a'jta^nec^  °nce  £  is  known. 

Now  for  x"^  we  have  from  convolution  arguments. 


n+1 

{  £  x„  P,  (i,n+l-j :H-j) 
j=l  J  1 


+  aix00  (i » n :M-1) j/p1  , 


(29) 


where  p^  is  the  probability  that  a  departing  job  is  interactive  and 


1 

x 

n 


/.  x 
i=0 


in 


M>n>0  . 


(30) 


Now  we  will  specify  the  condition  for  the  existence  of  Con¬ 
sidering  A,  it  is  clear  that  A  is  an  upper  triangular  matrix  whose 
entries,  with  notation  as  defined  in  the  previous  case,  are 


f  oo  M-i  "^l1"  i-i  "^i(M"jH 
of  (j-i)(1'e  e  dB0(t),  M>i>0;  M>j>i 


A(i,j)  =< 


0 


v. 


otherwise 


Clearly  A  is  reducible  and  therefore  the  condition  (4)  cannot  be 
directly  used.  However,  using  an  analysis  similar  to  that  of 
Lucantoni  [7]  page  7,  it  can  be  established  that  the  necessary 
and  sufficient  condition  for  stability  is 


xoy,po 


<  1  . 
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3.  NUMERICAL  RESULTS 


The  algorithm  generating  these  performance  measures  was  imple¬ 
mented  in  APL.  Some  sample  results  for  different  input  parameters 
Aq,  A^,  Uq,  y^,  and  M  for  exponential  service  times  are  shown  in 
Tables  1  and  2. 


Batch  Jobs 
Arrival 

Rate  AQ 

Mean  number  of  jobs  in  the  system  at  an  arbitrary 
time  epoch  when  higher  priority  is  given  to 

Interactive  Jobs 

Batch  Jobs  I 

Batch 

Interactive 

Batch 

Interactive 

.05 

0.2835 

0.5875 

0.2567 

0.6841 

.075 

0.5133 

0.7376 

0.4394 

0.9565 

.10 

0.8569 

0.8878 

0.6821 

1.2890 

.125 

1.4190 

1.0380 

1.0210 

1.6920 

.15 

2.4879 

1.1881 

1.5265 

2.1713 

.175 

5.2676 

1.3383 

2.3663 

2.7365 

Batch  jobs  service  rate  Pq  «  0.25 
Number  of  Interactive  terminals  N  **  5 
Interactive  job  arrival  rate  A^  *  0.1 
Interactive  jobs  service  rate  y^  *  2.0 


TABLE  1 


1 
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Number  of 
Interactive 
Terminals 

Mean  number  of  jobs  in  the  system  at  an  arbitrary 
time  epoch  when  higher  priority  is  given  to 

Interactive  Jobs 

Batch  Jobs  1 

Batch 

Interactive 

Batch 

Interactive 

5 

0.2835 

0.5875 

0.2567 

0.6841 

6 

0.2939 

0.7327 

0.2580 

0.8552 

7 

0.3038 

0.8895 

0.2591 

1.0407 

8 

0.3184 

1.0598 

0.2604 

1.2422 

Batch  jobs  arrival  rate  X^  =>  0.05 
Batch  jobs  service  rate  Uq  *  0.25 
Interactive  jobs  arrival  rate  X^  =  0.1 
Interactive  jobs  service  rate  *  2.0 


TABLE  2 

4.  CONCLUSION 

In  this  paper  we  have  modelled  the  computer  system  with  batch 
and  interactive  jobs  as  a  single  server  queueing  system  with  two 
classes  of  customers  and  non-preemptive  priority  service  disci¬ 
pline.  The  system  performance  measures  were  numerically  evaluated 
using  an  algorithmic  approach  proposed  by  Neuts  [9]  and  a  discrete 
state  level  crossing  analysis  (see  [11]). 
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APPENDIX 


(a)  Recursive  Evaluation  of  Transition  Probabilities 


Define 


P^(n,r:k)  *  P{n  batch  jobs  and  r  interactive  jobs  arrive  during 
a  time  interval  which  has  probability  distribution 
B^(*)  given  there  are  k  interactive  jobs  left  in 

the  source  at  the  beginning  of  this  time  interval} 


That  la 


PjCn.rik) 


o»e  °  (V>n  k  ~Xifc  r  "(k~r>xit 

Qf - ^ -  <p(l-e  1  )re  1  dB1(t),  i-0,1 


A0  k  »  “V  n  r  r  i  t 

(*></  e  t  (jr0(j)('1)Je  )e  dBl(t) 


X0  k  r  1  r  «  n  '<Xo+(k+J'r)Xlt 
rr  Q  t  (-Dj(I)  /  tne  0  1  dB.(t) 

nllj«0  3  0  1 


•  n 


(-1)°  nT  (r}  1  (-l>J(I)^iC)(xo+(k+j“r)xi)  ’ 


j-0 


where 


&5n*(a)  “  lim  ^.(s) 

1  s+a  dsn  1 


(Al) 
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and 


Vs)  =  0--dBi(x)  ,  i-0,1  . 
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Noce  that 


and 


P^n.rsk+l)  *  0  for  all  r>k+l 


Pi(n,0:k)  =  (-l)n  ^  &|n)  (A^kA^  ,  M>k>0  . 


Then  using  (A2),  P^(n,r:k)  can  be  calculated  in  a  recursive  fashion. 
That  la,  for  r«*l,  calculate  P^n.lsk),  k  *  1,2,...,M,  using  (A2)  and 
then  r»2,  etc. 


(b)  Quantities  Needed  In  the  Analysis  of  the  Imbedded  Markov  Chain 
(taken  from  Lucantoni  and  Neuts  [8]). 


A*(Z)  =  Z  ZrA 
r=0 


00 

JC(Z)  -  E  z\ 
r=*l 


A 

G(Z)  -  Z  ZA  Gr(Z)  =  ZA  (G(Z) ) 
r-0  r 


G  -  Z  A  Gr 


ire 

Wr 


1 
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! 


*  00  00  00  ~ 

JJ  -  (I-  Z  A  Gr_1)"1{e  +  [  I  a  -  E  A  Gr  l+  l  (r-l)A  GJd-frW)"1^} 

r-1  r  r-1  r  r-1  r  r-2  r  * 


00  oo 

{  m  ft  +  2  »  Gr“  Y  +  I  Z  V  *  BrGr_1  +  E  (r-l)BrG](I-0+^)’'1u 
r-1  r=l  r  r-1  r-2  r 


9 


MM 


1 


oo  00 

-lr  r-  „  „r-l 


-  jg  +  A0(I-B0)  \  +  {A0(I-B0)_'L[  E  B  -  E  B  G 

r-1  r-1 


OO  00 


+  E  (r-l)B  2]  +  E  A  -  E  A  Gr_1  +  E  (r-l)A  G  (I-&+G)"1* 
r-2  r  r-1  r  r-1  r  r-2  r  A 


A*n(Z)  *  -£  A*<Z> 

dz 


«(0)(1)  «  1,  4(0)<1)  -  X(0)(l)  -  TT,  S(1)  -  »A*(1)? 


jj(1)(i)  -  (i-a+h)"1!  -  5(1,s 


+ 

P 


§1 


i  4 
1  fi 


11 
i  a 
*  I 


?(1)W  “  ^(Dd-A+n)'1  -  6(1)£ 
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for  n>2 


l<«  -  -  "4MU)i(r), 


(n)  i  n 

jg  (i)  *  (i-A+n)-1  z 


(I-A+n)-1^(^)[A*(l)-6(r)(1)I]js(n-r)(1) 

'  ^WT)ah^a)k 


/nn>  a  T  /nuW/nr.* 


*  (1)  =  J „(r>rr'«H*n.t(i)  -  «(n-r)a)i!(i-*+n,-i 


(1)  1  “  « 

*  a> *  '  2^^  kf“h*  +  2t0  (»  *  to  j2Kk-l>V 

+  2«o  kf1kB^a>(1)  +  *o  k^y2,a>  - 

■  tlAot(  >(l>  +  ?(Ui5«<2)(1)}  - 

eV‘“1“ti0"  °  is  carried  °ut  “**“8  the  recursion  (see  page  6  of  [8]). 


G(0)  -  (I-Aj)"1^, 


G(k+1)  *  k>0 
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